Spin Transitions Induced by Magnetic Field in Quantum Dot Molecules 
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We present a theoretical study of magnetic field driven spin transitions of electrons in coupled 
lateral quantum dot molecules. A detailed numerical study of spin phases of artificial molecules 
composed of two laterally coupled quantum dots with N = 8 electrons is presented as a function 
of magnetic field, Zeeman energy, and the detuning using real space Hartree-Fock Configuration 
Interaction (HF-CI) technique. A microscopic picture of quantum Hall ferromagnetic phases cor- 
responding to zero and full spin polarization at filling factors v = 2 and v = 1, and ferrimagnetic 
phases resulting from coupling of the two dots, is presented. 
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I. INTRODUCTION 

The application of spin of electrons in quantum 
dots for generation of electron entanglement and quan- 
tum information processing in solid state devices is 
of current experimentali^^&i&^i&ii and theoretical 
interes^^^^ 6 -^^. Controlling the spin of electrons 
in single quantum dots by tuning the external magnetic 
field, the confining potential, number of electrons, and 
Zeeman coupling has been demonstrated&^ai&l. It was 
shown that in strong magnetic field electrons form a spin 
singlet quantum Hall droplet at filling factor v = 2. In- 
creasing the magnetic field leads to the spin-flip transi- 
tions until the spin polarized filling factor v = 1 droplet 
is reached. 2 Spin flips beyond the first spin flip are as- 
sociated with correlated states such as spin bi-excitons, 
identified and observed experimentally^.. Quantum dot 
molecules offer additional possibility of coupling and con- 
trolling spin transitions by tuning the tunneling barrier 
which controls the inter-dot coupling : 12 i 16 i 17 ' 18 i 19 The re- 
cently demonstrated time dependent control of the tun- 
neling barrier height and confining potential 6 , and the 
quantum state of the electron spin by applying oscillating 
magnetic field (Rabi oscillations )£, resulted in coherent 
manipulation of two electron spins in coupled quantum 
dot molecules. Recent experiments by Pioro-Ladriere et 
al. in RefH suggested that in strong magnetic field elec- 
trons are expected to form quantum Hall droplet in each 
quantum dot. Edge states of each droplet can be coupled 
in a controlled way using barrier electrodes, and at filling 
factor v = 2 effectively reduce the many-electron-double 
dot system to a two-level molecule^. When populated 
with one electron each, one expects to have singlet-triplet 
transitions of two valence electrons in the background of 
core electrons of the spin singlet v = 2 droplets. With 
increasing magnetic field, transitions to higher spin po- 
larized states are expected, with coupled quantum dots 
resembling artificial magnetic molecules. 

In this paper, we investigate the effect of the interdot 
tunneling and electron-electron interactions on the evolu- 
tion of total spin of electrons in a quantum dot molecule 
as a function of electron numbers and magnetic field. 



We study the many-body effects in the spin flip transi- 
tions by incorporating systematically the inter-dot and 
intra-dot electron-electron Coulomb interactions using 
real space Hartree-Fock Configuration Interaction (HF- 
CI) technique. We find quantum Hall droplets with zero 
and full spin polarization, identified as v — 2 and v = 1 
quantum Hall droplets^ 2 -, in analogy with single quan- 
tum dots and quantum Hall ferromagnets 20 . Between 
these two states, we find series of continuous transitions 
among partially spin polarized phases. These partially 
polarized phases correspond to spin flips. Simultaneous 
spin flip in each isolated dot must lead to even number 
of spin flips in a double dot. Recently, we have found 
partially spin polarized phases which correspond to odd 
number of spin flips^i in a double quantum dot. In Ref. 
pnj . we have identified these correlated states as quan- 
tum Hall ferrimagnets. 

Coherent superposition of two single particle levels in 
a double well potential in a form of symmetrical and an- 
tisymmetrical states is a well known example of quantum 
mechanical phenomena, from coherent charge oscillation 
of an electron between two states localized on two pro- 
tons in molecules, coherent oscillation between left 
handed and right handed Amino acids ,~ coherent control 
of Rabi oscillations of electron spin in quantum dots 7 , 
to macroscopic quantum resonance of Cooper pairs in 
mesoscopic superconducting grains 2 ^. The quantum Hall 
ferrimagnetic states, or spin unbalanced phases ; 2 1 25 are 
also a direct manifestation of coherent quantum mechani- 
cal tunneling and inter-dot electronic correlations. These 
states can be described in terms of linear combinations 
of spin excitons localized in left and right dots, which in 
turn lead to coherent spin oscillations, e.g., spin counter- 
part of coherent charge oscillations in molecules. 

The paper is organized as follows: In section [TT] we 
review the Hamiltonian of electrons confined in the lat- 
eral gated quantum dots. In section IIII1 the computa- 
tional methods of single particle configuration interaction 
(SP-CI) and unrestricted Hartree-Fock configuration in- 
teraction (URHF-CI) are summarized. To differentiate 
the spin transitions of quantum dot molecules and two 
isolated dots with zero inter-dot interaction, we briefly 
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present the spin phase diagram of single dots in sec- 
tion IIVI The microscopic picture of spin excitations in 
coupled quantum dots are discussed in sections [Vj and 
IVII The interpretation of spin excitations in terms of 
electron-hole excitations allows us to attribute the excita- 
tions with total spin S = 1 and S = 2 as spin exciton and 
bi-exciton. Pairing of excitons and the formation of bi- 
excitons due to strong inter-dot interaction is discussed. 
To elucidate the quantum Hall ferrimagnetic states, the 
real space representation of excitons is introduced. The 
paper is summarized in section I VIII 



II. HAMILTONIAN 

We describe electrons confined in quasi-two- 
dimensional quantum dots in a uniform perpendicular 
magnetic field by the effective mass Hamiltonian 

N 2 -I 

tf = £(T J + i^) + f-£-— (1) 



i=l 



where 



T=^(jV + - c A(f)) + V( X ,y) (2) 

is the single electron Hamiltonian in magnetic field. Here 
(f) = (x,y) describes electron position, A(r) = ^B x r 
is the vector potential, B is the external magnetic field, 
and V(r) is the quantum dots confining potential, m* 
is the conduction-electron effective mass, e is the elec- 
tron charge, and e is the host semiconductor dielectric 
constant (e = 12.8 in GaAs). Eiz = ^g/^B^izB is the 
Zeeman spin splitting, g is the host semiconductor g- 
factor (g = —0.44 in GaAs), /is is the Bohr magneton, 
and a is the Pauli matrix. In what follows, we present 
the numerical results in effective atomic unit (in GaAs 
effective Bohr radii = 9.79nm, and effective Rydberg 
Ry* = 5.93meV). 

The single particle eigenvalues (e^) and eigenvectors 
(ifi) are calculated by discretizing T in real space, and 
diagonalizing the resulting matrix using conjugate gradi- 
ent algorithms^. The details of this calculation can be 
found in Ref. H3. 



III. MANY BODY SPECTRUM 

To calculate the spectrum of interacting electrons, de- 
scribed by Hamiltonian H in Eq.|T]), we employ either 
the real space single particle or unrestricted Hartree- 
Fock states in configuration interaction techniques^ In 
the first SP-CI approach, single-particle levels arc used 
to construct many-electron configurations which are the 
basis of configuration interaction (CI) Hamiltonian. De- 
noting the creation (annihilation) operators for electrons 



in non-interacting SP state \aa) by S aa (c aa ), the Hamil- 
tonian of an interacting electron system in second quan- 
tization can be written as 

a a 

+ 2 ^ ] ^ 1 Vaa,/3a' j'ycr' ,/j,a c gcr c g a -i c 'y(r' c fj,(X (3) 
a/?7^i era' 

where the first term is the single par- 
ticle Hamiltonian, and V aa .f3 a > n<7 > ^ a = 
!drfdP(p%^cp^,(P)-^^(p^(?)cp^(f), is the 

two-body Coulomb matrix element. In the configuration 
interaction method the Hamiltonian of an interacting 
system is calculated in the basis of finite number of 
many-electron configurations. The total number of 
configurations (or Slater determinants participating in 
CI calculation) is determined by 



Nr. 



NJ 



NA(N S - 7V T ) 



NJ 



N^N.-Ny)] 



(4) 



Here N s is the number of single particle levels, and iVj 
and N± are the number of spin up and spin down elec- 
trons. This Hamiltonian is either diagonalized exactly 
for small systems or low energy eigenvalues and eigen- 
states are extracted approximately for very large number 
of configurations.— With increasing the number of single- 
particle levels N s , the number of configurations Nc grows 
very fast, yet a large number is needed to accurately ac- 
count for direct and exchange interaction, and electronic 
correlations. To improve the convergence of CI method 
we incorporate direct and exchange contribution into the 
basis states by replacing SP states with states obtained 
by the unrestricted Hartree-Fock method (URHF-CI). 
This implies expressing the new creation (annihilation) 
operators for URHF quasi-particles by a\ a (aw), with the 
index i representing the URHF orbit quantum numbers. 
The URHF basis can be expanded in a linear combina- 
tion of SP states. In terms of SP creation (annihilation) 
operators we write 



N, 



(5) 



(i) 

The transformation coefficients, X aa , satisfy the self- 
consistent Pople-Nesbet equations 



18,28,29,31. 



£{ £ A,+ £aw 7 ££A;^Ag, 

7=1 a, /3=1 a' j=l 

-wEE^w}^ = € F a«, (6) 

a' j=l 

where e^ F are the URHF eigenenergies. The TV-lowest 
energy URHF levels form a Slater determinant occupied 
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by HF quasi-electrons, corresponding to HF ground state. 
The rest of orbitals with higher energies are outside of 
the HF Slater determinant (unoccupied HF levels), con- 
tribute to electronic correlations, and can be used for CI 
calculation. The many body Hamiltonian of the inter- 
acting system in the URHF basis can finally be written 
as: 



H 



^2^{ip\ T \jp) a \* a i<? 

a ij 

ijkl u a' 

where ^Ji a ,j a ' ,ha' ,ltr are the Coulomb matrix elements in 
the URHF basis' Here 



(ia\T\ja) 



(ia\U H - Ux\j(r), 



(8) 



where Uh and Ux are the Hartree and exchange opera- 
tors 

(ia\U H \ja)= £ ^ 7 ,A;WA«EE A -' )A ?^ 

a(3fj,iy— 1 a' k—\ 



N, 



N„ 



ia\U x \ja)= ]T K^A^Ag^A^AW. 

a/3 [iv — 1 



k=l 



The resulting CI Hamiltonian matrix constructed in the 
basis of URHF configurations is either diagonalized ex- 
actly for small systems, or low energy eigenvalues and 
eigenstates are extracted approximately for very large 
number of configurations^!. The details of the calcu- 
lation and the convergence of the results as a function 
of number of basis and configurations, and comparison 
between SP-CI and URHF-CI methods can be found in 
Ref. M 
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FIG. 1: (Color online) Single particle spectrum as a func- 
tion of cyclotron frequency uj c for a gaussian single dot with 
strength Vo = — lORy* , and A = 2.5ao in the presence of 
Zeeman splitting. Arrows represent spin of electrons. For il- 
lustration purposes a very high Zeeman coupling g = —9 is 
used. 
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FIG. 2: (Color online) The energy of spin configurations 
shown in the boxes, \v — 2), \X), and \XX) with total spin 
S = 0, S = 1, and S = 2, using LLL orbitals, and Ez = 
0. E v= 2 (S — 0) is the reference of energy. The arrows 
surrounded by circles represent the holes. 



IV. SPIN TRANSITIONS IN A SINGLE DOT 



We describe a single dot by an isotropic gaussian con- 
fining potential V(x,y) = Voe~( x +v ''^ . The single 
particle eigenenergies of such quantum dot as a function 
of cyclotron frequency uj c = eB/m* are shown in Fig. [T] 

States with spin up (down) are shown by bold (dashed) 
lines. For illustration, a very large g-factor is introduced. 
We note in Fig. [TJ with increasing magnetic field the en- 
ergy of spin up (bold) levels decrease. These levels, and 
their spin down partners, correspond to the levels of the 
lowest Landau level (LLL). We now populate the lowest 
energy states with a number of electrons. From previous 
work^, the minimum number of electrons which exhibits 
all nontrivial phenomena in the spin evolution of a single 
quantum dot is N — 4. The N = 4 configurations which 
minimize the kinetic and Zeeman energy are shown in 
Figffl Due to the crossing of spin up and down levels, 



there are three different configurations S z = 0, S z = 1, 
and S z = 2. These configurations illustrate increasing 
spin polarization of the electronic droplet with increas- 
ing magnetic field. With very small Zeeman energy the 
increasing spin polarization in quantum dots is driven by 
electron-electron interactions. We hence turn off the Zee- 
man coupling and turn on electron-electron interactions. 

We start with the lowest energy configuration build 
with SP LLL states, the S z — spin singlet v = 2 con- 
figuration \v = 2). The spin excitations with S = 1 and 
S = 2, are constructed by removing electrons from oc- 
cupied states and putting to unoccupied states, and can 
be described in terms of single exciton \X) and biexciton 
\XX) configurations, as shown in Fig. [21 Neglecting the 
mixing between configurations, we calculate the energy of 
each spin configuration. The result of this calculation is 
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FIG. 3: (Color online) The spin evolution of the ground state 
of single quantum dot as a function of magnetic field using 
URHF-CI method with N 3 = 8 HF levels, corresponding to 
Nc = 784 configurations. 

shown in Fig. O We chose the energy of the 5 = 0, v = 2 
state {E v= 2) as the reference energy. As shown in Fig. 
with increasing magnetic field even without Zeeman 
energy both the first and second spin flip transitions oc- 
cur at uj c — 1.6 and lo c — 2.75 due to electron-electron 
interactions. 

The effect of correlations on the evolution of spin of 
electrons in a single quantum dot obtained by the URHF- 
CI method is shown in Fig. [3] Here the spin unpolar- 
ized (5 Z — 0) URHF states have been constructed out 
of Ni = 30 SP states. N s = 8 HF levels have be taken 
to construct CI Hamiltonian, resulting in Nc = 784 con- 
figurations. At zero magnetic field we find 5 = 1 triplet 
due to Hund's rule for electrons in a half-filled p-shell. 
With increasing magnetic field, the single particle energy 
gap opens up, leading to suppression of the 5 = 1 triplet 
state, and formation of the v = 2 (5 = 0) singlet state. 
The first spin flip 5 = 1 state appears around co c = 3, fol- 
lowed by the second spin flip spin polarized 5 = 2 state 
at to c — 4.5. The flipping of the second spin is interrupted 
by a low-spin, 5 = 0, strongly correlated state. This state 
was previously identified with the formation of spin sin- 
glet bi-exciton£ The first and second spin flip state can 
be obtained both for the noninteracting electrons and in 
Hartree-Fock approximation while the formation of spin 
singlet bi-exciton is a result of electronic correlations in 
a quantum dot. 



V. SPIN TRANSITIONS IN QUANTUM DOT 
MOLECULES 

We now turn to study the spin transitions of laterally 
coupled quantum dots. We describe the molecule by elec- 
tron (Nl,Nr) and ground state spin numbers (Sl,Sr) 
of isolated left (L) and right (R) dots. The spin phase 
diagram turns out to depend on electron numbers in each 
dot. Here we will focus on molecules build out of identical 



dots with Nl = Nr. For a given number of electrons the 
magnetic field will be used to tune their individual spin 
Sl = Sr. The goal will be to determine the total spin 
of the molecular system. The molecular coupling will be 
controlled by the height of the tunneling barrier. To il- 
lustrate the physics we will discuss in detail quantum dot 
molecule (4, 4) with four electrons each and contrast its 
properties with a single N = 4 quantum dot discussed in 
previous section. 



A. One electron spectrum of a quantum dot 
molecule 

We parameterize quantum dot molecule po- 
tential in terms of a sum of three Gaussians 
V(x,y) = V L eM~ (x+a L W ] + Vr exp[- + 
V p exp[— 4 ■ Here Vl,Vr describe the depth 

P- Py 

of the left and right quantum dot minima located at 
x = —a,y = and x = +a, y — 0, and V p is the 
central plunger gate potential controlling the tunneling 
barrier. The confining potential is parameterized as 
V L = V R = -10, A = 2.5, a = 2, and A Px = 0.3, 
Ap y = 2.5, in effective atomic units. The single particle 
eigenvalues and eigenfunctions are calculated numeri- 
cally by discretization of the Schrodinger equation with 
the quantum dot molecule potential. The parameters 
of the confining potential considered in this work are 
chosen to represent weakly coupled quantum dots. 
At zero magnetic field the SP levels of the electrons 
exhibit well separated S,P, and D electronic shells, and 
at high magnetic field they form molecular shells of 
closely spaced pairs of bonding- antibonding orbitals. 
The half-filled molecular shells correspond to electron 
numbers (N L = 2k - 1, N R = 2k - 1) and filled shells 
correspond to (Nl — 2k, Nr — 2k) configurations (k 
is integer). The resulting single particle spectrum as a 
function of magnetic field has been presented recently 
in Fig. 1 of Ref. |2l|. In order to understand and 
visualize the spectrum in high magnetic field we expand 
confining potential in the vicinity of each minimum. 
This gives a parabolic potential V(r) = m* io^t 1 j2 with 
the strength ojq — 2y/\Va\/A 2 . The low energy spectrum 
of each dot corresponds to two harmonic oscillators with 
eigen-energies e nm = uj + (n + 1/2) + w_(m + 1/2). Here 
lu± = yWg + /4 i u) c /2, lo c is the cyclotron energy, 
and n, m = 0, 1, 2, .... With increasing magnetic field the 
lo- decreases to zero while a>+ approaches the cyclotron 
energy lu c , and the states \m, n) evolve into the nth 
Landau level. In high magnetic field the corresponding 
wavefunctions admit a description in terms of localized 
LLL orbitals 18 . In this limit linear combinations of the 
LLL orbitals m from left and right dot forms molecular 
shells of closely spaced symmetric-antisymmetric pairs 
with eigen-energies expressed approximately as 
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e m \ a =(jJ-[m+-)- A— 



(9) 



Here Ez — 7^ c is the Zeeman energy, a = +1 (—1) corre- 
sponds to spin | (|), 7 = ™*g and A is the pseudospin in- 
dex: the symmetric (antisymmetric) orbitals are labeled 
by A = +1 (—1), the parity of the molecular orbitals. 
A m is the symmetric-antisymmetric gap. The Zeeman 
coupling induces spin splitting with increasing magnetic 
field. This has been illustrated in Fig. 1 of Ref. [U We 
now populate the quantum dot molecule with N = 8 
electrons. This is an example of electronic configurations 
corresponding to filled molecular shells. Fig. 1 of Ref. I2H 
shows the evolution of the lowest energy states of nonin- 
teracting electrons, with the corresponding total spin 5 
states separated by vertical lines. We find S = 0, 2, and 
4 phases with even 5 which correspond to simultaneous 
spin flips in each isolated dot. 

However, we also find odd 5 phases. The first odd 
spin state with 5 = 1 occurs between magnetic fields 
corresponding to uo* cl m 3.25 and lo* 2 rs 3.9 where 



=2,ct=T.A=+1 



= l,,7 = l,X = 



and 



=2,ct=T,A= 



e m =i, CT =i,A=+i- The odd spin flip is related to the split- 
ting of energy levels due to tunneling. Using single par- 
ticle eigen-energies given by Eq.®, we find the first spin 
flip at the value of the magnetic field B\ 



7W C (1) = w-(l) 



A 2 (l) + A 1 (l) 



(10) 



where Zeeman energy equals the single dot level splitting 
minus the average symmetric- antisymmetric gap for the 
two levels involved. The second spin flip takes place at 
w c {2) such that 



7U/ C (2) = w_(2) + 



A 2 (2) + Ai(2) 



(11) 



Hence the difference in the magnetic fields corresponding 
to first and second spin flips is a direct measure of the 
tunneling splitting: 



7(w c (2) - w c (l)) 



A 2 (2) + Aj(2) 
2 

, A 2 (l) + Ai(l) 



(12) 



From the spectrum of Fig. 1 of Ref. [2l|, we observe that 
the states with odd spins 5=1, and 5 = 3, are sta- 
ble within narrow range of magnetic fields due to spin 
flip transitions among the electrons that occupy the lev- 
els with energy separation proportional to the inter-dot 
tunneling amplitude. This is in contrast with the first 
spin flip transition in single dots (compare with Fig. [1]) 
which is stable in a wide range of magnetic fields. For this 
reason the existence of odd spin states in the spin phase 
diagram of quantum dot molecules can be interpreted as 
the measure of inter-dot interaction. 



B. Many electron quantum dot molecule spectrum 

In this section we present an analysis of spin transitions 
driven by electron-electron interaction. We focus on the 
tunnel coupled lowest Landau level orbitals m. Denoting 
the creation (annihilation) operators for electron in non- 
interacting SP state \m\(j) by o| nAcr (c m \ a ) (with a as 
spin label), the Hamiltonian of an interacting system, 
Eq. [TJ can be written as 



H 



EE 



tmXaCjnXfj Cm Act 



2 D D' 

{m.A} 00' 



m\\\<j, 1712X2^' \ V \m3X3a' , m^Xicr) 



X C m 1 A 1 a C m 2 A 2 a' C " l 3 A 3 <r' c «u A 4 cr 



(13) 



The single particle states of coupled quantum dot 
molecules in magnetic field are labeled by the orbital 
quantum numbers m, and the pseudospin index A. The 
first term in Eq. [13] is the single particle Hamiltonian, 
and V a /3(iv is the two-body Coulomb matrix element. 
Here {a, /?, n, v} represent the states with quantum num- 
bers (m, A, a). We now turn to the construction of the 
relevant configurations. 



C. S = 0, v = 2 state 

The v = 2 state of quantum dot molecule with TV elec- 
trons and total spin 5 = 0, shown in Fig. 2J is the 
product of electron creation operators 

N/4-1 

i"= 2 >= n n n c Lj«> m 

m=0 A=l,2er=T,l 

The energy associated with this state 
N/4-1 

m=0 A=±l <t 

can be expressed in terms of self-energy S(m, A, a) 

N/4-1 

£(m, A, a) = E E ( 2 ( mX > m'X'\V\m'X', mX) 

m'=0 A'=±l 

-{mX,m'X'\V\mX,m'X')). (16) 



D. 



5=1 spin exciton 



The 5=1 spin flip excitation is constructed by remov- 
ing an electron from occupied j orbital in a v = 2 state, 
and putting it into an unoccupied state i, 



2), 



(17) 
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FIG. 4: (Color online) The basis of spin configurations in high 
magnetic fields. The first spin transition states S = 1 identify 
with two independent set: {Xss, Xaa}, and {Xsa, Xas}- In 
the former the electron-hole transitions occurs between the 
states with the same symmetry and hence they do not mix 
with the latter which exhibit the process of electron-hole ex- 
citations between states with opposite parity. 



where j = (to, A, J,), and i = (to', A', | ). Denoting quasi- 
particle energy levels (electrons dressed by interaction) 
by Si = £j + the energy of one exciton follows 



AEx^, = £i -£j - (i,j\V\j,i) 



(18) 



where ABxj^ = Ex,^ — E v —2 is the energy of exciton 
relative to the = 2 state energy. The last term in Eq. 
(|18p is the electron-hole Coulomb attraction. The lowest 
energy states of the single exciton corresponding to the 
first spin flip state are depicted in Fig. 2J We classify 
the single excitons based on their relative parity with re- 
spect to the parity of the v = 2 state. Pair of states 
(Xss, Xaa) describes transitions between pairs of levels 
with the same parity and so parity is conserved by these 
transitions. In contrast, spin flip transitions represented 
by (Xsa, Xas) do not conserve parity as they describe 
transitions between pairs of levels with opposite parity. 
We identify (Xss, Xaa), and (Xsa, Xas) by their parity 
quantum numbers n = +1 and n = — 1 , respectively. It 
is important to note that Coulomb interactions do not 
mix excitons with different parities, and the CI Hamil- 
tonian constructed in the basis of spin excitons is block 
diagonal. These pairs of independent excitons are shown 
inside the boxes in Fig. [U 

In real space, single excitons can be expressed in terms 
of linear combination of excitons localized in each dot. 
This basis is shown in Fig. \5\ In each isolated dot, at 
a critical field, a transition from the S = singlet to 
the 5 = 1 triplet state takes place. This configuration, 
corresponding to Xrr and Xll in Fig. is equivalent 
to a localized electron-hole excitation in one dot in the 
presence of background of electrons in the other dot. Be- 
cause of the geometrical symmetry associated with the 
electron-hole excitation, the ground state of the system 
without external bias and interdot tunneling has double 
degeneracy: The state with (Sl = 0,Sr = 1) has exactly 
the same energy as the state (Sl = 1, Sr — 0). The 
many body wavefunction of such a molecular state can 
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FIG. 5: (Color online) The lowest energy excitons in real 
space \X RR ) = |X( lii j a) ^(2 ifl , T )>, \X LL ) = \X {1>l ,d^(2,l,-\)), 
\X RL ) = |X(i i ij ji )_ > ( 2 ,£,t)), and \X LR ) = \X^ l ^^ 2 ,r,-[)) ■ 
The magnetic ordering of these states can be described by 
means of ferrimagnetic coupling (Sl = 0,Sr — 1) for Xrr, 
and (Sl = 1,Sr = 0) for Xll, and ferromagnetic coupling 
(Sl = 1/2, Sr = 1/2) for X RL , and X LR . 



be expressed as linear combination of degenerate states 
\S L = 0,S R = 1} and \S L = 1,S R = 0). We identify 
these pairs of molecular states with quantum Hall ferri- 
magnets. For a range of magnetic field these two states 
are separated from another pair of single excitations Xlr 
and X RL with (Sl = 1/2,Sr = 1/2) by an energy gap 
due to Coulomb interactions. In general, in the case of 
filled shells, molecular states with odd total spin S cor- 
respond to quantum Hall ferrimagnets. Our analysis of 
spin transitions in real space is presented in section IVT1 



E. S = 2 spin bi-exciton 

With increasing magnetic fields, a higher polarized 
state with S = 2, (Sl = 1, Sr = 1) equivalent to a spin 
bi-exciton state appears. A bi-exciton is constructed by 
removing a pair of electrons from occupied states and 
putting them into unoccupied states (see Fig. 0]), 



\XXj^ k ^i) = c\c]cjC k \u = 2) 



(19) 



The energy of biexciton can be decomposed into the en- 
ergy of two single excitons plus their interaction 



AE 



AE 



A£^, + 6V, 



(20) 



where AEx-^i has been defined in Eq. [T5] SV is the 
binding energy between two excitons, accounting for the 
electron-electron, electron-hole and hole-hole interactions 



SV = (l,i\V\i,l) - (l,i\V\l,i) 

-(l,j\V\j,l)-(k,i\V\i,k) 
+(j,k\V\k,j)-(j,k\V\j,k). 



(21) 
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FIG. 6: (Color online) (a) The energies of two single spin 
excitons Xsa and Xas with odd parity, the energy of the 
odd parity correlated exciton Xsa+as, and the energy of the 
spin bi-exciton S = 2 state as a function of magnetic field. 
All energies measured from the energy of the v = 2, S = 
state, (b) The energies of two single spin excitons Xss and 
Xaa with even parity, the energy of the even parity correlated 
exciton Xss+aa, and the energy of the spin bi-exciton 5 = 2 
state as a function of magnetic field. All energies measured 
from the energy of the v = 2, S = state. 



F. First versus second spin flip in a quantum dot 
molecule 

Unlike in a single quantum dot, first spin flip state 
5 = 1 is not an eigenstate of the coupled quantum 
dot Hamiltonian. There are two possible spin excitons 
for a given parity, and they are coupled by Coulomb 
interactions. We use two distinct single exciton basis 
{X S s, X AA ] and {X S a, X as }, labeled by parity n = +1 
and 7r = — 1, to construct the two 2x2 Hamiltonians. We 
denote by AE Xss +x AA , AE XsA +x AS the lowest eigen- 
energies of these Hamiltonians. 

Fig. (Ha) shows the numerically calculated energies of 
odd parity spin excitons as a function of magnetic field. 



The energy AEx AS of spin exciton Xas is positive for 
magnetic fields shown but the energy AEx SA of spin ex- 
citon Xsa becomes negative at uj c = 3.8 i.e. the Xsa 
spin flip state becomes the lower energy state than the 
v = 2, 5 = state. However, in stark contrast with a 
single quantum dot, Fig. [21 we find that the second spin 
flip state XX becomes the ground state at lower mag- 
netic field lj c = 2.1. Hence, unlike in a single quantum 
dot we find a transition from spin singlet 5 = state 
directly to 5 — 2 second spin flip state. This is a transi- 
tion corresponding to even total spin numbers, as if the 
two dots were flipping their spin simultaneously. How- 
ever, correlations or mixing of the two spin excitons Xas 
and Xsa lowers the energy of the spin exciton. The en- 
ergy AEx SA +x AS ) °f the correlated single spin state is 
significantly lower and equals both the energy AExx of 
the second spin flip bi-exciton and of the v = 2, S = 
state at cu c — 2.1. At this value of the magnetic field the 
energy of the bi-exciton and of the exciton are almost 
identical, and we might expect that for larger number of 
configurations correlations will stabilize the single spin 
flip exciton even further. The effect of correlations on 
the even parity excitons Xaa and Xss is shown in Fig. 
EJb) . We see that mixing of the two even parity excitons 
lowers the energy AEx ss +x AA of the correlated single 
spin state. This energy equals both the energy AExx of 
the second spin flip bi-exciton and of the v — 2, 5 = 
state at uj c = 2.1. By comparison with Fig. [6{a) we see 
that the value of co c = 2.1 also corresponds to the change 
in parity of the single spin exciton 5 = 1 state. 

As illustrated in Fig. [51 with mixing of single exciton 
configurations the energy of quantum dot molecule ex- 
hibits four-fold degeneracy at where 5 = 0, 5 = 2 
and two different parity 5=1 states become the low- 
est energy states. The 5=1 states show stability in 
a narrow range of magnetic field, within the accuracy 
of our numerical results. With further increase of mag- 
netic field, single excitons condense into pairs of excitons 
forming bi-excitons. The existence of single, odd, spin 
excitons is hence a signature of electronic correlations. 
These states do not exist at the Hartree-Fock level. 

To support the assertion that correlations are respon- 
sible for the existence of odd spin excitons we employ 
URHF-CI calculation. From the solution of the Pople- 
Nesbet equations we obtain HF eigen-energies, shown in 
Fig. [7| (a), for TV = 8 electrons with S z = 0, as a function 
of the magnetic field. The HF wavefunctions are used as 
a basis in HF-CI calculation of the ground state. We em- 
ploy N s — 8 HF basis states (equivalent to Nc = 4900 
configurations) to calculate the total spin of electrons in 
quantum dot molecule as a function of the magnetic field. 
From this calculation we find that the 5 = 0, 5 = 1, and 
5 = 2 states are almost degenerate at u)* x = 3.1. The 
prediction of URHF-CI is in qualitative agreement with 
effective SP-CI model, presented above. The direct, ex- 
change, and correlation energies calculated by URHF-CI 
shift the transition point to higher magnetic fields. 

With increasing of the magnetic field we find that the 
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FIG. 8: (Color online) (a) The evolution of spin of the N = 8 
electron quantum dot molecules as a function of magnetic field 
for g = and V p = 7 (solid line). For comparison the spin 
evolution of two isolated dots is shown (dashed line). The 
width of odd spin plateaux have been artificially enlarged to 
be visible by eyes. The vertical dashed line (with purple- 
color online) tends to show qualitatively a range of u c in 
which the gap vanishes. The horizontal dashed line (with 
red-color online) shows the corresponding spin one state with 
zero gap. (b) The evolution of the energy gap of the quantum 
dot molecules as a function of magnetic field. 



FIG. 7: (Color online) (a) URHF energies vs. cyclotron en- 
ergy for N = 8 electrons with S z = 0. (b) Evolution of lowest 
energies for S = 0, 1, 2 states calculated by URHF-CI method. 
The energies are measured from the energy of the S = state. 



gap between different total spin states tends to vanish. 
This gapless phase is seen in Figs. J7| and © in the 
vicinity of lo c = 4. This phase is followed by the 5 = 2, 
5 = 3 and 5 = 4 states. The latter state corresponds to 
a fully spin polarized v = 1 droplet. 

Let us now summarize the similarities and the differ- 
ences in the evolution of total spin of two isolated dots 
and a quantum dot molecule. Fig. [8] (a) shows the evo- 
lution of total spin with increasing magnetic field for two 
nonintcracting N — 4 quantum dots and for N = 8 quan- 
tum dot molecule. Fig[5fb) shows the energy gap of the 
molecule as a function of magnetic field. The noninter- 
acting quantum dots spin evolution is obtained by adding 
results from two isolated dots, each dot evolving with 
magnetic field according to Fig. [31 We find that the ef- 
fect of interdot interaction is to renormalize the magnetic 
fields at which spin transitions take place, and more im- 
portantly, to lead to the appearance of odd spin states 
5 = 1 and 5 = 3. While the existence of odd spin states 
is most striking, the presence of spin polarized phases 
is also nontrivial. The fact that spins of electrons on 
two quantum dots align demonstrates the existence of 
ferromagnetic dot-dot coupling. In the case of antifer- 




FIG. 9: (Color online) The effect of increasing Zeeman energy 
on the (a) evolution of spin and (b) evolution of energy gap 
of the N = 8 electron quantum dot molecules as a function 
of magnetic field. The width of g = odd spin plateaux have 
been artificially enlarged to be visible by eyes. 



romagnetic coupling there would have been no net spin 
eventhough each dot has finite spin. We will show in the 
next section that such antiferromagnetic coupling exists 
in the N = 8 molecule at low magnetic fields. 

Finally, following the gap in FigJSJb) we find the ex- 
istence of low spin "gapless" states, molecular analogs 
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FIG. 10: (Color online) The effect of detuning, i.e., difference 
in confining potential between left dot Vl and right dot Vr 
on the (a) evolution of spin and (b) evolution of energy gap 
of the N = 8 electron quantum dot molecules as a function 
of magnetic field. 



of spin singlet bi-excitons. Hence Fig[5] shows that elec- 
tronic inter-dot correlations stabilize the odd spin phases 
but their stability range is very narrow. As shown in Sec. 
(IV Ap without electron-electron interactions competition 
between quantum mechanical tunneling and Zeeman en- 
ergy was responsible for the existence of odd spin phases. 
The effect of finite Zeeman energy is similar in an inter- 
acting system. FiglH] shows the effect of increasing Zee- 
man energy on the evolution of spin and energy gap of the 
N = 8 electron quantum dot molecules as a function of 
magnetic field. All parameters are the same as in FigJS] 
We see that increasing Zeeman energy renormalizes the 
magnetic field value of spin flips and, more importantly, 
stabilizes the odd spin phases. From the evolution of the 
energy gap shown in Fig ^jp we also see the vanishing of 
the low spin depolarized phase in the vicinity of lu c = 4 
and the stabilization of the spin polarized S — 2 phase. 
The Zeeman energy depends on the g- factor. For GaAs 
the g-factor is gcaAs — —0.44 while for InAs and InSb the 
g-factors are g InAs = -14, and g InSb = -BoM^L^s 
Hence by adding In one can hope to tune the g-factor of 
quantum dot molecules. 

To conclude our analysis of the N — 8 electron quan- 
tum dot molecule we discuss the effect of asymmetry be- 
tween the two dots. While for molecules built out of two 
atoms each component is identical, quantum dots are de- 
fined by gates or etching and one must understand the 
effect of differences between the two dots on the stability 
diagram 30 . In Fig[llj]we show the evolution of spin and 
energy gap of the N — 8 electron quantum dot molecules 
as a function of magnetic field. The two dots are differ- 
ent, with confining potential of the left dot Vl = —10 
unchanged but the potential of the right dot detuned by 
IRy to Vr = —11. As anticipated, the effect of detun- 



ing results in increased stability of the odd spin flip state 
5 = 1. 



VI. REAL SPACE ANALYSIS OF SPIN 
TRANSITIONS IN QUANTUM DOT 
MOLECULES 

In this section, we describe the analysis of spin transi- 
tions in real space. Eq. (flU)) describes the Hamiltonian 
of an interacting system in second quantization in non- 
interacting SP state | to A<r). The single particle states 
of coupled quantum dot molecules in magnetic field are 
labeled by the single dot orbital quantum numbers to, 
and the pseudospin index A. The symmetric (antisym- 
metric) orbitals are labeled by A = 1 (—1), the parity 
of the orbitals in two symmetric dots, to represents 
the combined Landau level, and angular momentum 
quantum numbers, to = (n, I). The first term in Eq. (|13p 
is the single particle Hamiltonian, and V aa ^ a i ^ va = 
J dfj dr'^ aa (r)^ a , (P)-^.^ (r')^(f), is the 

two-body Coulomb matrix element. Here {a,/3, /i, v\ 
represent the states with orbital quantum numbers 
(to, A). 

Alternatively denoting the creation (annihilation) op- 
erators for electron in non-interacting localized SP state 
\mscr) by <i mso . (d msa ), the Hamiltonian of an interacting 
system in second quantization can be written as 



ms a 

+ ^ t t m ^ ] (1 — 3sis 2 )d mSl<J d, 1 



ma S1S2 



- ^2 y^raisi 0- . ■m 2 s 2 cr / \V\m 3 s 3 a', m 4 s 4 a) 



{m,s} <t<t' 

-ft 



Yfft d 1 d ,d 

A "mi si tr u m 2 s 2 a' u "m 3 s 3 a' «m 4 s 4 er 



(22) 



Here s = 1 (2) are pseudospin labels of electron localized 
in left (right) dot. The relation between Eqs. HH and 
| can be established by a rotation in pseudospin space 



J = i 

-mA<7 ^2 



rTsE^iA- 1 ^"- We find e ™ = (Cm - A = 1 



£m,A=-i)/2, t m — {e m ,\=i — £m,A=-i)/2, and 

(m\S\a, m 2 S2a'\V\m 3 s 3 a' , TO 4 s 4 cr) = - ^ A^ 1 " 3 ^ 2 " 1 

A 

\ s 3 3 - 1 \ s i l - 1 (m 1 \ 1 (T,m 2 \2<7'\V\m 3 \ 3 a',m i \ 4 <7) (23) 

A. 5 = ground state 

The v = 2 state of quantum dot molecule with N elec- 
trons and total spin S = is the product of spin polarized 
localized electrons 



N/4,-1 

= 2 )= n n n 

m=0 s=l,2a=1,l 



(24) 
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The energy associated with this state follows 

AT/4-1 2 

E u=2 = E ( 2g "» + s ( m . *)) (25) 

m— s—1 

where E(m, s) is the electron self-energy: 

N/4-1 2 

S(m, s) = (2(ms, m! s\V\m s , ms) 



m'— s' — 1 



(ms, m' s'\V\ms, m's'}) 



(26) 



B. S = 1 exciton 

In each isolated dot, at a critical field, and driven 
by electron-electron Coulomb interaction and increasing 
electron kinetic energy, a transition from S — singlet to 
S = 1 triplet is seen. The cost in kinetic energy is lowered 
if localized electrons in one dot flip the spin (Sl = 1), 
while the other electrons in the second dot occupy the 
lowest energy single particle states to form spin singlet 
droplet (Sr = 0). This configuration corresponding to 
Xrr (or Xll) in Fig. [5l is equivalent to a localized 
electron-hole excitation in one dot in the presence of the 
background of electrons in the other dot. Because of the 
geometrical symmetry associated with the electron-hole 
excitation, the ground state of the system (without any 
external bias and interdot tunneling) has double degen- 
eracy: The state with (Sl = 0,Sr = 1) has exactly 
the same energy as the state (Sl — ^,Sr — 0). The 
many body wavefunction of such a molecular state can 
be expressed as linear combination of degenerate states 
\Sl = 0,S r = 1), and \S L = 1,S R = 0). We identify 
these pairs of excitations as quantum Hall fcrrimagnets. 
For a range of magnetic fields these two states are sep- 
arated from another pair of single excitations Xlr and 
X RL with (S L = 1/2, Sr = 1/2) by an energy gap due 
to Coulomb interaction. 

An excitonic state is constructed by removing an elec- 
tron from occupied state and putting into an unoccupied 
state 



\X j ^ i )=d\d j \v = 2) 



(27) 



where j = (m, s,J.), and i = (m',s',|) The lowest en- 
ergy basis of the single exciton (first spin flip state) is 
depicted in Fig. [5] Labeling the direct and indirect 
spin flip transitions (with ferrimagnetic and ferromag- 
netic spin ordering) by {Xrr,X L l}, and {X LRl X RL } 
and using their symmetries we find that the direct (in- 
direct) states are two-fold degenerate Ex RR = Ex LL 
(Ex LR = Ex RL )- Here the subscripts are defined as 
RR = (l,R,i) -> (2,i?,|), LL = (l,L,i) -» (2,L,|), 
RL = (1,R, i) -f (2, L, j), and LR = (1, L, |) -> (2, R, | 
). Note that in non- interacting system the basis is four 
fold degenerate E Xrr = E Xll = E Xlr = E Xrl - 



Denoting quasi-particle energy levels by 
the energy of one exciton reads 



E 



E u 



(i,j\V\j,i) 



e,- + S(0 



(28) 



where the last term is the electron-hole Coulomb inter- 
action. In the basis of single exciton states, the Hamil- 
tonian of the QD molecules can be expressed as 

reff 



rrcff T-reff 

-"4x4 — J 4x4 



4x4; 



where 



J 4x4 




(29) 



(30) 



is the non-interacting part, and 



V A 



cff 

4x4 



Vrrll Vrrlr Vrrrl 
Vf?Rr.r. 11 VllLR VllRL 

Vlrrl 




RR.LL 

V* V* 

v RR.LR. v LLLR 

V* V* V* 



(31) 



is the Coulomb interaction matrix between the single ex- 
citon states. 

Note that the pair of states Xrr and Xll is not cou- 
pled by the single tunneling term because the scattering 
process between Xrr and Xll requires the exchange of 
two particles simultaneously. For that reason this pro- 
cess is second order in tunneling. The same is true for 
the states Xlr, and Xrl- 



C. 5 = 2 bi-exciton 

With increasing magnetic fields, a higher polarized 
state with 5 = 2, (Sl = 1,Sr = 1) equivalent to a 
bi-excitonic state tends to appear as ground-state. A bi- 
exciton is constructed by removing a pair of electrons 
from occupied states and putting into unoccupied states 



IX,- 



.,k->i) 



2). 



(32) 



The energy of biexcitonic state can be decomposed into 
the energy of two single excitons plus their interaction 



AE 



= AE 



AE 



Xk—<-i 



sv. 



(33) 



SV is the binding energy between two excitons, accounted 
for the electron-electron, electron-hole and hole-hole in- 
teractions 

SV = (Z,*|V|t,0 - (l,t|V|/,<) 
~{l,3\V\j,l) 
-(k,i\V\i,k) 

+{j,k\V\k,j)-{j,k\V\j,k) (34) 

In the quantum dot molecule considered in this study, we 
find SV < for large magnetic fields, i.e., two isolated 
excitons favor to pair and form a biexcitonic state where 
the energy of a biexciton is lower than energy of two 
isolated single excitons. 
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FIG. 11: The energies of spin excitons with respect to v — 
2, S — state as a function of magnetic field. Ef s is the 
lowest eigen-energy of the effective Hamiltonian iS|x4- The 
inset illustrates the energy levels of single particle localized 
states. 



D. Excitonic condensation: a SP-CI effective model 



energy of biexciton (5 — 2), calculated from the energy 
of the ground state with 5 = 0, are shown in Fig. [Sj 

As predicted by single configuration SP-CI, the di- 
rect spin flip transition takes place at frtu*i ~ 2.8 where 
E Xrr < E v —2, as shown in Fig. [TTJ Within this range 
of magnetic fields, the energy of a biexciton is lower than 
the energy of a single exciton due to strong Coulomb in- 
teraction. At hu* 2 ~ 2.1, a transition to S = 2 state 
is seen due to pairing of single spin excitons. However, 
because of strong mixing between the single excitonic 
states, electron correlations improves the energy of the 
5 = 1 state, enough to bring the first spin flip transition 
point to fiu>*i « fi.w*2, where three states with 5 = 0, 
5 = 1, and 5 = 2 appeared to be almost degenerate. 
The Zeeman coupling remove such degeneracy. As a re- 
sult, the 5 = 1 state tends to become stable in a narrow 
range of magnetic field. 

VII. CONCLUSION 



The eigenvalues of T|^4 follow E AS < E ss < E AA < 
Esa ( a s t 2 > ti). The corresponding eigenstates are 
\X AS ),\X SS ), \Xaa), I Asa), shown in Fig. H Note that 
for non-interacting electrons AE X = E x — E v —2 > 
where E° x = Ex RR = E Xll = E Xrl = E Xlr . By 
simplifying the Coulomb interaction among electrons as 
Vijki — > Vq > 0, we can calculate the self-energy, and 
the ground state energy of 5 = state analytically as: 
£ = NVo/2, and E v = 2 = E° =2 + (N/2) 2 V . The latter 
is the energy of single exciton measured from the ground 
state energy. We can also calculate the energy of sin- 
gle exciton AE x . ( = — ej — Vq < AE X , and the 
energy of the biexciton AE XX = 2AE X — 2Vq, where 
AE X = AE Xrr = AE Xll , and SV = -2V < 0. The 
energy difference between biexciton and a single exciton 
follows AE XX - AE X = AE X - 2Vo. In the limit of 
strong Coulomb interaction (large magnetic fields) we 
find ti — 6j < 3Vb and E xx < E x . 

The energies of two distinct excitons corresponding to 
the direct and indirect first spin flip transition, Xrr, 
and Xlr, the lowest energy eigenvalue of Hf^ A , and the 



In conclusion, we have presented detailed analysis of 
the magnetic field driven spin transitions in quantum 
dot artificial molecules with N = 8 electrons as a func- 
tion of external magnetic field, Zeeman energy, and the 
detuning, using Hartree-Fock configuration interaction 
method. The magnetic field allows the tuning of the to- 
tal spin of electrons in each artificial atom. Quantum 
mechanical tunneling and electron-electron interactions 
couple spins of each artificial atom and result in ferro- 
magnetic, anti-ferromagnetic , and ferrimagnetic states 
of quantum dot artificial molecules, tunable by the mag- 
netic field and barrier potential. 
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